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Abstract 

We point out that resonance saturation in QCD can be understood in the 
large- N c limit from the mathematical theory of Pade Approximants to meromor- 
phic functions. These approximants are rational functions which encompass any 
saturation with a finite number of resonances as a particular example, explaining 
several results which have appeared in the literature. We review the main prop- 
erties of Pade Approximants with the help of a toy model for the {VV — AA) 
two-point correlator, paying particular attention to the relationship among the 
Chiral Expansion, the Operator Product Expansion and the resonance spectrum. 
In passing, we also comment on an old proposal made by Migdal in 1977 which 
has recently attracted much attention in the context of AdS/QCD models. Fi- 
nally, we apply the simplest Pade Approximant to the (VV — AA) correlator 
in the real case of QCD. The general conclusion is that a rational approximant 
may reliably describe a Green's function in the Euclidean, but the same is not 
true in the Minkowski regime due to the appearance of unphysical poles and/or 
residues. 



1 Introduction 



The strong Chiral Lagrangian is a systematic organization of the physics in powers of 
momenta and quark masses, but requires knowledge of the low-energy constants (LEC) 
to make reliable phenomenological predictions. As with any other effective field theory, 
these LECs play the role of coupling constants and contain the information which comes 
from the integration of the heavy degrees of freedom not explicitly present in the Chiral 
Lagrangian (e.g. meson resonances). 

At 0(p 4 ) there are 10 of these constants [lj. Although at this order there is enough 
independent information to extract the values for these constants from experiment, 
this will hardly ever be possible at the next order, 0(p 6 ), because the number of 
constants becomes more than a hundred [5]. In the electroweak sector the proliferation 
of constants appears already at 0(p 4 ) \3\. Although in principle these low-energy 
constants may be computed on the lattice, in practice this has only been accomplished 
in a few cases for the strong Chiral Lagrangian at 0(p A ), and only recently [I]. 

The large N c expansion [5] stands out as a very promising analytic approach capa- 
ble of dealing with the complexities of nonperturbative QCD while, at the same time, 
offering a relatively simple and manageable description of the physics. For instance, 
mesons are qq states with no width, the OZI rule is exact and there is even a proof 
of spontaneous chiral symmetry breaking [6]. Furthermore, its interest has recently 
received a renewed boost indirectly through the connection of some highly supersym- 
metric gauge theories to gravity [7], although the real relevance of this connection for 
QCD still remains to be seen. However, in spite of all this, the fact that no solution to 
large- N c QCD has been found keeps posing a serious limitation to doing phenomenol- 
ogy. For instance, in order to reproduce the parton model logarithm which is present 
in QCD Green's functions in perturbation theory, an infinity of resonances is necessary 
whose masses and decay constants are in principle unknown. 

On the other hand, QCD Green's functions seem to be approximately saturated by 
just a few resonances; a property which has a long-standing phenomenological support 
going all the way back to vector meson dominance ideas [8], although it has never 
been properly understood. In a modern incarnation, this fact translated into the very 
successful observation [5] that the strong LECs at 0(p 4 ) seem to be well saturated by 
the lowest meson in the relevant channelp after certain constraints are imposed on some 
amplitudes at high-energy in order to match the expected behavior in QCD [TTJ [12]. 
It was then realized that all these successful results could be encompassed at once as 
an approximation to large- N c QCD consisting in keeping only a finite (as opposed to 
the original infinite) set of resonances in Green's functions. This approximation to 
large- N c QCD has been termed Minimal Hadronic Approximation (MHA) [T3] because 
it implements the minimal constraints which are necessary to secure the leading non- 
trivial behavior at large energy of certain Green's functions through the marriage of 
the old resonance saturation and the large- iV c approximation of QCD. In recent years, 
a large amount of work has been dedicated to studying the consequences of these ideas 

m. 

1 This is less clear in the scalar channel, however. See Ref. [TO] 
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However, the high-energy matching with a finite set of resonances, first suggested in 
[TT] , makes it clear that the treatment is not amenable to the methods of a conventional 
effective field theory. An effective field theory is an approximation for energies smaller 
than a heavy particle's mass and, therefore, cannot deal with momentum expansions 
at infinity as in the case of the Operator product Expansion (OPE). In other words, 
the fact that the set of resonances in each channel is really infinite precludes the 
naive expansion at large momentum because there is always a mass in the spectrum 
which is even larger. The sum over an infinite set of resonances and the expansion 
for large momentum are operations which do not commute [15]. In those Green's 
functions containing a contribution from the parton model logarithm, this is made 
self-evident since a naive expansion at large momentum can only produce powers and 
not a logarithm, which is why large- iV c QCD requires an infinity of resonances in the 
first place. 

The problem can be delayed one power of a s if one requires the use of the resonance 
Lagrangian [9] to be limited only to Green's functions which are order parameters of 
spontaneous chiral symmetry breaking. These order parameters vanish to all orders in 
a s in the chiral limito and, therefore, avoid the presence of the parton model logarithm 
which, otherwise, would preclude from the outset any matching to a finite number of 
resonances. However, the concept of a Lagrangian whose validity is restricted only to a 
certain class of Green's functions has never been totally clear; and even if the resonance 
Lagrangian is restricted by definition to order parameters, the problem surfaces again 
in the presence of logarithmic corrections from nontrivial anomalous dimensions, which 
make the exact matching at infinite momentum impossible. 

In a slightly different context, a somewhat similar observation was also made in Ref. 
|16j . In this paper it was observed that it is impossible to satisfy the large momentum 
fall-off expected in large- N c QCD for the form factors which can be defined through 
a three-point Green's function, if the sum over resonances in the Green's function is 
restricted to a finite set. Interestingly, this again pointed to an incompatibility of the 
QCD short- distance behavior with an approximation to large N c which only kept a 
finite number of resonances. 

A further piece of interesting evidence results from the comparison between the 
analysis in Refs. [17] and [TT] . After imposing some good high-energy behavior in 
several Green's functions and form factors including, in particular, the axial form 
factor governing the decay n — > jeU, Ref. [UJ obtains, keeping only one vector state V 
and one axial-vector state A, that their two masses must be constrained by the relation 
Ma = \/2My. The work in Ref. [TT], on the contrary, does not use the axial form factor 
and obtains, after performing a very good fit within the same set of approximations, 
the precise values My = 775.9 ± 0.5 MeV and M A = 938.7 ± 1.4 MeV. These values 
for the masses, although close, are not entirely compatible with the previous relation. 
In other words, the short- distance constraint from the axial form factor is not fully 
compatible with the short-distance constraints used in [T7J if restricted to only one 
vector and one axial-vector states^. 

2 E.g., the two-point correlator (VV — AA). 

3 Adding one further state does not change the conclusion |17j . 
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In this paper we would like to suggest that all the above properties can be un- 
derstood if the approximation to large N c QCD with a finite number of resonances is 
reinterpreted within the mathematical Theory of Pade Approximants (PA) to mero- 
morphic functions |18j . QCD Green's functions in the large N c and chiral limits have 
an analytic structure in the complex momentum plane which consists of an infinity 
of isolated poles but no cut, i.e. they become meromorphic functions [19]. As such, 
they have a well-defined series expansion in powers of momentum around the origin 
with a finite radius of convergence given by the first resonance masfl This is all that 
is needed to construct a Pade Approximant. A theorem by Pommerenke [20] assures 
then convergence of any near diagonal PA to the true function for any finite momen- 
tum, over the whole complex plane, except perhaps in a zero-area set. The poles of 
the original Green's function (i.e. the resonance masses) belong to this zero-area set 
because not even the original function is defined there, but there are also extra poles. 
These extra poles are called "defects" in the mathematical literature [18]. When the 
Green's function being approximated is of the Stieltjes type^l, the poles of the PA are 
always real and located on the Minkowski region Re(g 2 ) = Re(— Q 2 ) > 0, approaching 
the physical poles as the order of the PA is increased [2T]. However, this takes place 
in a hierarchical way and, while the poles in the PA which are closest to the origin are 
also very close to the physical masses, the agreement quickly deteriorates and one may 
find that the last poles are several times bigger than their physical counterparts [22]. 
The same is true of the residues. In section 3, we will see with the help of a model 
that the same properties are met in a meromorphic function whose spectral function 
is not positive definite, except that some of the poles in the PA may even be complex. 

This means that Minkowskian properties, such as masses and decay constants, 
cannot be reliably determined from a PA except, perhaps, from the first poles which 
are closest to the origin. If not all the residues and/or masses are physical, then there 
is no reason why they should be the same in the form factor governing 7r — > ^ev and 
in the Green's function (VV — AA), explaining the different results found in [TT] and 
[TT] we alluded to above. Furthermore, the form factors of all but the lightest mesons, 
defined through the residues of the corresponding 3-point Green's functions, will not 
be reliably determined from a PA to that Green's function, again in agreement with 
the findings in [16J. 




The situation in the Euclidean is different. In general, PAs cannot be expanded at 
infinite momentum to generate an OPE type expansion for the true function. Never- 
theless, Pommerenke's theorem assures a good approximation at any finite momentum, 
no matter how large. Of course, the order of the PA will have to increase, the larger 
the momentum region one wishes to approximate. For instance, in Ref. [21] it was 
shown with the help of a simple model how, even in the case of the (VV) correlator 
which contains a logQ 2 at large values of Q 2 > 0, the PAs are capable of approximat- 
ing the true function at any arbitrarily large (but finite) value of Q 2 > 0, without the 
need for a perturbative continuum. In section 3 we will show, again with the help of a 

4 The pion pole can always be eliminated multiplying by enough powers of momentum. We are 
assuming here the existence of a nonvanishing gap in large- N c QCD. 

5 Roughly this means that the associated spectral function is positive definite, like in the case of 
the two-point correlator (VV). See Ref. [18] (chapter 5) for a more precise definition. 




3 



model, how this is also true in the more general case of a non-positive definite spectral 
function such as (VV — AA). This means that PAs are a reliable way to approximate 
the original Green's function in the Euclidean but not in the Minkowskian regime. 

In 1977, A. A. Migdal [23] suggested PAs as a method to extract the spectrum of 
large-iV c QCD from the leading term in the OPE of the (VV) correlator, i.e. from 
the parton model logarithm. However, nowadays this proposal should be considered 
unsatisfactory for a number of reasons [21], the most simple of them being that different 
spectra may lead to the same parton model logarithm [25J. In fact, the full OPE series 
is expected to be only an asymptotic expansion at Q 2 = oo (i.e. with zero radius 
of convergence), and PAs constructed from this type of expansions cannot in general 
reproduce the position of the physical poles [26]. For instance, we show this explicitly 
with the help of a model for (VV — AA) in the Appendix. Migdal's approach has 
been recently adopted (in disguise) in some models exploiting the so-called AdS/QCD 
correspondence [27] and, consequently, the same criticism also applies to them. 

In Ref. [28J a model for the (VV — A A) two-point correlator with a spectrum 
consisting of an infinity of resonances was suggested as a theoretical laboratory for 
studying the relationship between the spectrum and the coefficients of the OPE. In 
this paper several conventional methods usually employed in the literature were tested 
against the exact result from the model. These included: Finite Energy Sum Rules as 
in Ref. [33], pinched weights as in Ref. [31], Laplace transforms as in Ref. |35j and, 
finally, also resonance saturation as in the MHA method. The bottom line was that 
no method was able to produce very accurate predictions for the OPE coefficients. In 
all the methods but the last one, the reason for this lack of accuracy was basically due 
to the fact that the OPE requires an integral over the whole spectrum, whereas the 
integral is actually cut off at an upper limit (in the real case, the upper limit is m T ). 
This is why even if one uses the real spectrum the result may be inaccurate [29]. In 
the case of the MHA the reason was, as we will comment upon below, that the poles 
were not allowed to be complex. 

In section 3 we will revisit this (VV — AA) model, now from the point of view of 
PAs. The model reproduces the power behavior of QCD at large Q 2 > except that 
the model is simple enough not to have any log Q 2 and, therefore, it cannot reproduce 
the nonvanishing anomalous dimensions which exist in QCD. We do not think this is 
a major drawback, however, because in QCD these logarithms are always screened by 
at least one power of a s and, hence, in an approximate sense, it may be licit to ignore 
them. In the model such an approximation becomes exact^. Will the PAs be able to 
reproduce the large Q 2 expansion of the (VV— AA) model? We will see that the answer 
is affirmative. Therefore, the reason why the MHA method was not able to predict 
accurately the OPE coefficients in Ref. [28] is because the lowest PA has complex 
poles which were not allowed in |28) . When these complex poles are considered, the 
accuracy achieved is better and, most importantly, improves for a higher PA. Since 
the model allows the construction of PAs of a very high order, we have checked this 
convergence up to the Pade , which is able to reproduce the first non vanishing 
coefficient of the OPE in the model with an accuracy of 52 decimal figures. Together 

6 For a model with a logQ 2 , the reader may consult Ref. [21]. 
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with other numerical examples which will be discussed in section 3, we take this as a 
clear evidence of the convergence of the method. This renders some confidence that 
PAs may also do a good job in the real case of QCD. 

The rest of the paper is organized as follows. In section 2 we review some generali- 
ties of rational approximants, in section 3 we describe the (VV — AA) model and apply 
different rational approximants to learn about the possible advantages and disadvan- 
tages of them. In section 4 we apply the simplest PA to the case of the real {VV — AA) 
two-point function in QCD. Finally, we close with some conclusions. 



2 Rational approximations: generalities 

Let a function f(z) have an expansion around the origin of the complex plane of the 
form 

oo 

/(*) = , (!) 

n=0 

One defines a Pade Approximant (PA) to f(z) , denoted by P^f(z), as a ratio of two 
polynomials Qm{z), Rn{z)j, of order M and N (respectively) in the variable z, with 
a contact of order M + N with the expansion of f(z) around z = 0. This means 
that, when expanding Ppf(z) around z = 0, one reproduces exactly the first M + N 
coefficients of the expansion for f(z) in Eq. flTJ: 

P%{*) = ^44 « fo + h z + h z 2 + ... + W z M+N + 0(^ +M+1 ) • (2) 

R N {z) 

At finite z, the rational function P^f \z) constitutes a resummation of the series (00). Of 
special interest for us will be the case when N = M + k, for a fixed k, because then the 
function behaves like l/z k at z = oo. The corresponding PAs P^j +k {z) belong to what 
is called the near-diagonal sequence for k ^ 0, with the case k = being the diagonal 
sequence. 

The convergence properties of the PAs to a given function are much more difficult 
than those of normal power series and this is an active field of research in Applied 
Mathematics. In particular, those which concern meromorphic functional are rather 
well-known and will be of particular interest for this work. The main result which 
we will use is Pommerenke's Theorem [20J which asserts that the sequence of (near) 
diagonal PA's to a meromorphic function is convergent everywhere in any compact set 
of the complex plane except, perhaps, in a set of zero area. This set obviously includes 
the set of poles where the original function f(z) is clearly ill-defined but there may 
be some other extraneous poles as well. For a given compact region in the complex 
plane, the previous theorem of convergence requires that, either these extraneous poles 
move very far away from the region as the order of the Pade increases, or they pair 
up with a close- by zero becoming what is called a defect in the mathematical jargon 



7 Without loss of generality we define, as it is usually done, Rn(0) = 1. 

8 A function is said to be meromorphic when its singularities are only isolated poles. 



5 



[30] . These are to be considered artifacts of the approximation. Near the location of 
these extraneous poles the PA approximation clearly breaks down but, away from these 
poles, the approximation is safe. 

In the physical case the original function f(z) will be a Green's function G(Q 2 ) 
of the momentum variable Q 2 . In QCD in the large N c limit this Green's function 
is meromorphic with all its poles located on the negative real axis in the complex Q 2 
plane. These poles are identified with the meson masses. On the other hand, the region 
to be approximated by the PAs will be that of euclidean values for the momentum, i.e. 
Q 2 > 0. The expansion of G{Q 2 ) for Q 2 large and positive coincides with the Operator 
Product Expansion. 

In general a meromorphic function does not obey any positivity constraints and, 
as we will see, this has as a consequence that some of the poles and residues of the 
PAs may become compleo^. This clearly precludes any possibility that these poles and 
residues may have anything to do with the physical meson masses and decay constants. 
However, and this is very important to realize, this does not spoil the validity of the 
rational approximation provided the poles, complex or not, are not in the region of 
Q 2 one is interested in. It is to be considered rather as the price to pay for using a 
rational function, which has only a finite number of poles, as an approximation to a 
meromorphic function with an infinite set of poles. 

When the position of the poles in the original Green's function is known, at least 
for the lowest lying states, it is interesting to devise a rational approximation which has 
this information already built in. The corresponding approximants are called Partial 
Pade Approximants (PPAs) in the mathematical literature [31 J and are given by a 
rational function ¥^ K (Q 2 ): 

pAf (n 2\ _ Qm(Q 2 ) /„s 

nAQ } ~ Rn(Q 2 ) t k (Q 2 ) ' {6) 

where Qm(Q 2 ), Rn(Q 2 ) and Tk{Q 2 ) are polynomials of order M, N and K (respec- 
tively) in the variable Q 2 . The polynomial Tk(Q 2 ) is defined by having K zeros 
precisely at the location of the lowest lying poles of the original Green's functional 

T K (Q 2 ) = (Q 2 + Ml) (Q 2 + Ml) ... (Q 2 + Ml) . (4) 

As before the polynomial Rn{Q 2 ) is chosen so that i2jv(0) = 1 and, together with 
Qm{Q 2 ), they are defined so that the ratio k(Q 2 ) matches exactly the first M + N 
terms in the expansion of the original function around Q 2 = 0, i.e. : 

Kk(Q 2 ) ^h + hQ 2 + h Q 4 + ■■■ + fu + N Q 2M+2N + o{Q 2N+2M+2 ) . (5) 

At infinity, the PPA in Eq. © obviously falls off like 1/Q 2N+2K ~ 2M . Exactly as 
it happens in the case of PAs, also the PPAs will have complex poles for a general 



A special case which does obey positivity constraints is when the function is Stieltjes. In this 
case the poles and residues of the PAs are purely real and with the same sign as those of the original 
function [2T] , 

10 For simplicity, we will assume that all the poles are simple. 
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meromorphic function, which prevents it from any interpretation in terms of meson 
states. 

Finally, another rational approximant defined in mathematics is the so-called Pade 
Type Approximant (PTA) [31J T^(Q 2 ) : 

npMf^2\ _ Qm(Q 2 ) / k \ 

where T^(Q 2 ) is also given by the polynomial (jlj), now with N preassigned zeros at 
the corresponding position of the poles of the original Green's function, G(Q 2 ). The 
polynomial Qm{Q 2 ) is defined so that the expansion of the PTA around Q 2 = agrees 
with that of the original function up to and including terms of order M + 1, i.e. 

T^(Q 2 ) ^fo + fiQ 2 + h Q 4 + - + fu Q 2M + 0(Q 2M+2 ) . (7) 

At large values of Q 2 , one has that TjJ(Q 2 ) falls off like 1/Q 2N - 2M . Clearly the PTAs 
particular case of the PPAs, i.e. T^(Q 2 ) = F*f N (Q 2 ) and coincide with what has 
been called the Hadronic Approximation to large- N c QCD in the literature |13j . 

Let us summarize the mathematical jargon. A Pade Type Approximant (PTA) 
is a rational function with all the poles chosen in advance precisely at the physical 
masses. A Pade Approximant (PA) is when all the poles are left free. The intermediate 
situation, with some poles fixed at the physical masses and some left free, corresponds 
to what is called a Partial Pade Approximant (PPA). 



3 Testing rational approximations: a model 

Let us consider the two-point functions of vector and axial-vector currents in the chiral 
limit 

n^(g) = i J d 4 xe^{J^ A (x)4 V > A (0)) = (q,q v - g, v q 2 ) H^(g 2 ) , (8) 

with Jy(x) = d(x) , ~f tJ- u(x) and J^(x) = d(x)'~f tl '~f 5 u(x). As it is known, the difference 
nv(<? 2 ) ~ ^A(q 2 ) satisfies an unsubtracted dispersion relation^] 

LV_ A (g 2 ) = / ^— - ImLV_ A (t) . (9) 

Jo t - q 2 - ie 7T 

Following Refs. [32l [28], we define our model by giving the spectrum as 

1 oo 

-ImLV(t) = 2F 2 5(t-M 2 )+2y^F 2 (n)5(t-M 2 (n)) , 

n=0 

1 oo 

-Imn A (t) = 2F 2 5(t) +2V F 2 (n)5(t - M 2 A (n)) . (10) 

7T L — ' 

n=0 

11 The upper cutoff which is needed to render the dispersive integrals mathematically well defined 
can be sent to infinity provided it respects chiral symmetry [15] . 



7 



Here F p , M p are the electromagnetic decay constant and mass of the p meson and 
FyA(n) are the electromagnetic decay constants of the n — th resonance in the vector 
(resp. axial) channels, while My^(n) are the corresponding masses. F is the pion 
decay constant in the chiral limit. The dependence on the resonance excitation number 
n is the following: 

Fy A (n) = F 2 = constant , M VA (n) = m VA + n A 2 , (11) 

in accord with known properties of the large-iV c limit of QCD [5] as well as alleged 
properties of the associated Regge theory [37J. 
The combination 

n LR (q 2 )= 1 -(Tl v ( q 2 )-Tl A ( q 2 )) (12) 

thus reads 

2 F[ F 2 ~ f F 2 F 2 

LR[q ' ~ q 2 + -q 2 + M 2 + ^ [ -g 2 + M 2 (n) -g 2 + M\ (n) 

This two-point function can be expressed in terms of the Digamma function ip(z) = 
£\ogV(z) as [28] 

n , 2S F 2 F 2 F 2 f/-g 2 + ^i\ ,f-q 2 +ml\} 

n ^^# + 3^+^{n^^J"H^^J}- (14) 

To resemble the case of QCD, we will demand that the usual parton-model logarithm 
is reproduced in both vector and axial-vector channels and that the difference (Q has 
an operator product expansion which starts at dimension six. A set of parameters 
satisfying these conditions is given b}@ 

F = 85.8 MeV , F p = 133.884 MeV , F = 143.758 MeV , (15) 
M p = 0.767 GeV, m A = 1.182 GeV, m v = 1.49371 GeV , A = 1.2774 GeV , 

and is the one we will use in this section. This set of parameters has been chosen to 
resemble those of the real world, while keeping the model at a manageable level. For 
instance, the values of F p and M p in (fT5]) are chosen so that the function Ulr in (IT41) 
has vanishing 1/Q 2 and 1/Q A in the OPE at large Q 2 > 0, as in real QCD. In fact, the 
model admits the introduction of finite widths (which is a 1/N C effect) in the manner 
described in Ref. [32], after which the spectral function looks reasonably similar to 
the experimental spectral function. This comparison can be found in Fig. 5 of Ref. 
[28] . But this model is also interesting for a very different reason. In Ref. [28] several 
attempts were made at determining the coefficients of the OPE by using the methods 
which have become common practice in the literature. Among those we may list Finite 
Energy Sum Rules [33], with pinched weights [34] , Laplace sum rules [35J and Minimal 
Hadronic Approximation [13]. As it turned out, when these methods were tested on 

12 These numbers have been rounded off for the purpose of presentation. Some of the exercises which 
will follow require much more precision than the one shown here. 
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Co 


c 2 


c 4 


c 6 




C_ 6 


C-8 


-7.362 


21.01 


-43.92 


81.81 


-2.592 


1.674 


-0.577 



Table 1: Values of the coefficients C 2 k from the high- and low-Q 2 expansions of 
Q 2 TI L r(-Q 2 ) in Eq. (TJ2P in units of 10~ 3 GeV 2 ~ 2k . Notice that C_ 2 = and 
Co = —Fq (the pion decay constant in the chiral limit), see text. 



the model, none of them was able to produce very accurate results. We think that this 
makes the model very interesting (and challenging !) as a way to assess systematic 
errors |36j . 

Defining the expansion of the Green's function ([9]) in Q 2 = — q 2 around Q 2 = 0, oo 

as 

Q 2 U LR (-Q 2 ) « C 2 k Q 2k , with fc = 0,±l,±2,±3,... (16) 



one obtains that the coefficients accompanying inverse powers of momentum, akin to 



the Operator Product Expansion at large Q 2 
k — 1 — p) : 



> 0, are given by (p = 1,2,3,... with 



C 



2k 



-F 2 S pA + 



NP+l 



F 2 M 2p ~ 



P 



B,, 



my 
A 2 



— B r , 



A 2 



(17) 



As stated above, F p and M p are defined 



where B p (x) are the Bernoulli polynomials [TO 
by the condition that the above expression ( fTTj) vanishes for k = 0, — 1 enforcing that 
Q 2 Hlr(—Q 2 ) ~ Q~ 4 at large momentum, as in QCD. We emphasize that the above 
coefficients of the OPE in Eq. (1171) can not be calculated by a naive expansion at large 
Q 2 of the Green's function in Eq. fflBl . In other words, physical masses and decay 
constants do not satisfy the Weinberg sum rules [15]. 

On the other hand, for the coefficients accompanying nonnegative powers of mo- 
mentum, akin to the chiral expansion at small Q 2 , one has (k = 1, 2, 3, ...): 



Co = ~F ( 



C 



2k 



:-i) 



fc+i 



M 2k 



(k - 1)! A 2k 



(k-i) 



my 
A 2 



A 2 



where il)^ k ~ l \z) = d k ~ 1 ip(z) /dz k ~ 1 . In Tabled] we collect the values for the first few of 
these coefficients C 2 k- 

Let us start with the construction of the rational approximants to the function 
Q 2 Ulr{—Q 2 )- Since our original function (fill) falls off at large Q 2 as Q~ A , this is a 
constraint we will impose on all our approximants. 

The simplest PA satisfying the right falloff at large momentum is P^iQ 2 ), so we 
will begin with this case. In order to simplify the results, and unless explicitly stated 
otherwise, we will assume that dimensionful quantities are expressed in units of GeV 
to the appropriate power. Fixing the three unknowns with the first three coefficients 
from the chiral expansion of (1141) (i.e. Co 2 4) one gets the following rational function 



R 



(Q 2 + zr)(Q 2 + z* R ) 



R 



3.379 xl0~ 3 , z R = 0.6550 + i 0.1732 . (19) 



9 



Im(q 2 ) 



n 

10 1 H 




Figure 1: Location of the poles (dots) and zeros (squares) of the Pade Approximant P^i—Q 2 ) 
in the complex q 2 plane. We recall that Q 2 = —q 2 . Notice how zeros and poles approximately 
coincide in the region which is farthest away from the origin. When the order of the Pade 
is increased, the overall shape of the figure does not change but the two branches of complex 
poles move towards the right, i.e. away from the origin. 



We can hardly overemphasize the striking appearance of a pair of complex-conjugate 
poles on the Minkowski side of the complex Q 2 plane. Obviously, this means that these 
poles cannot be interpreted in any way as the meson states appearing in the physical 
spectrum ( TT0]fT3T) . In spite of this, if one expands ( |T9l) for large values of Q 2 > 0, one 
finds C_4 = — r\ = —3.379 x 10~ 3 which is not such a bad approximation for this 
coefficient of the OPE, see Table HJ Even better is the prediction of the fourth term in 
the chiral expansion, which is Cq = 79.58 x 10~ 3 . 

This agreement is not a numerical coincidence and the approximation can be sys- 
tematically improved if more terms of the chiral expansion are known. In order to 
exemplify this, we have amused ourselves by constructing the high-order PA P^iQ 2 )- 
This rational approximant correctly determines the values for C_4 5 _6,-8 with (respec- 
tively) 52, 48 and 45 decimal figures. In the case of C103, which is the first predictable 
term from the chiral expansion for this Pade, the accuracy reaches some staggering 192 
decimal figures. This is all in agreement with Pommerenke's theorem [20] . 

As it happens for the PA (TT91) . also higher-order PAs may develop some artificial 
poles. In particular, Figure [D shows the location of the 52 poles of the PA P^iQ 2 ) m 
the complex q 2 plane. Of these, the first 25 are purely real and the rest are complex- 
conjugate pairs. A detailed numerical analysis reveals that the poles and residues 
reproduce very well the value of the meson masses and decay constants for the lowest 
part of the physical spectrum of the model given in d 1 3ti 151) . but the agreement deteri- 
orates very quickly as one gets farther away from the origin, eventually becoming the 
complex numbers seen in Fig. [IJ It is by creating these analytic defects that rational 
functions can effectively mimic with a finite number of poles the infinite tower of poles 
present in the original function (fT4"j) . 

For instance the values of the first pole and residue in P^iQ 2 ) reproduce those 
of the p in ffTBT) within 193 astonishing decimal places for both. However, in the case 
of the 25th pole, which is the last one still purely real, its location agrees with the 
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physical mass only with 3 decimal figures. This is not to be considered as a success, 
however, because after the previous accuracy, this is quite a dramatic drop. In fact, 
the residue associated with this 25th pole comes out to be 29 times the true value. 
The lesson we would like to draw from this exercise should be clear: the determination 
of decay constants and masses extracted as the residues and poles of a PA deteriorate 
very quickly as one moves away from the origin. There is no reason why the last 
poles and residues in the PA are to be anywhere near their physical counterparts and 
their identification with the particle's mass and decay constant should be considered 
unreliable. Clearly, this particularly affects low-order PAs. 

A very good accuracy can also be obtained in the determination of global euclidean 
observables such as integrals of the Green's function over the interval < Q 2 < oo. 
Notice that the region where one approximates the true function is far away from the 
artificial poles in the PA. For instance, one may consider the value for the integral 

POO 

h = (-1) / dQ 2 Q 2 U LR (Q 2 ) = 4.78719 x 10~ 3 , (20) 
Jo 

which, up to a constant, would yield the electromagetic pion mass difference in the 
chiral limit [38] in the model (fl4j) . The PA P^iQ 2 ) reproduces the value for this 
integral with more than 42 decimal figures. This suggests that one may use the integral 
( 1201 as a further input to construct a PA. 

For example if we fix the three unknowns in the PA P 2 °(Q 2 ) by matching the first 
two terms from the chiral expansion but now we complete it with the pion mass dif- 
ference (l20j) instead of a third term from the chiral expansion as we did in (JTSJ the 
approximant results to be 

^i.Q 2 ) = (n2 , i m > with r R = 2 ' 898 x 10_3 > Z R = °- 5618 + 1 °' 2795 ' 

(Q 2 + z R ){Q 2 + Z* R ) 

(21) 

This determines C_ 4 = —2.898 x 10 3 and C 4 = —41.26 x 10 3 , which shows that using 
the pion mass difference is not a bad idea. Notice how the position of the artificial pole 
has changed with respect to f[T^]) . 

Artificial poles and analytic defects are transient in nature, i.e. they appear and 
disappear from a point in the complex plane when the order of the Pade is changed. On 
the contrary, the typical sign that a pole in a Pade is associated with a truly physical 
pole is its stability under these changes in the order of the Pade. Of course, when the 
order in the Pade increases there have to be new poles by definition, and it is natural 
to expect that some of them will be defects. Pade Approximants place some effective 
poles and residues in the complex Q 2 plane in order to mimic the behavior of the true 
Green's function, but it can mimic the function only away from the poles, e.g. in the 
Euclidean region. Obviously, PAs cannot converge at the poles, in agreement with 
Pommerenke's theorem [20J, since not even the true function is well defined there. The 
point is that what may look like a small correction in the Euclidean region may turn 
out to be a large number in the Minknowski region. To exemplify this in simple terms, 

13 We remark that this procedure, although reasonable from the phenomenological point of view, 
strictly speaking lies outside the standard mathematical theory of rational approximants [HI [31] . 
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let us consider a very small parameter e and imagine that a given Pade P(Q 2 ) produces 
the rational approximant to the true Green's function G(Q 2 ) given by 



G{Q 2 ) « P{Q 2 ) = R{Q 2 ) + > ( 22 ) 

where R{Q 2 ) is the part of the Pade which is independent of e. Although for Q 2 > 
there is a sense in which the last term is a small correction precisely because of the 
smallness of e, for Q 2 < this is no longer true because of the pole at Q 2 = —M 2 . 
This pole is in general a defect and may not represent any physical mass. In fact, 
associated with this pole, there is a very close-by zero of the Pade P{Q 2 ) at Q 2 = 
—M 2 — e R(—M 2 ) , as can be immediately checked in (1221) . This is another way of 
saying that a defect is characterized by having an abnormally small residue and is the 
origin of the pairs of zeros and poles in the y-shaped branches of Fig. HJ Therefore, 
not only are defects unavoidable but one could say they are even necessary for a Pade 
Approximant to approximate a meromorphic function with an infinite set of poles. 

Similarly to masses, also decay constants may be unreliable. To see this, imagine 
now that our Pade is given by 

p e 

P ( Q ) = Q2 + M 2 + (Q2 + M 2) (Q2 + M 2 + e 2 ) ' ^ 

again for a very small e. As before, the term proportional to e may be considered 
a small correction for Q 2 > 0. However, at the pole Q 2 = —M 2 the decay constant 
becomes F + e _1 which, for e small, may represent a huge correction. When the poles 
are preassigned at the physical masses, like in the case of PTAs, it is the value of the 
residues that compensates for the fact that the rational approximant lacks the infinite 
tower of resonances. As we saw before, the residues of the poles in the Pade which lie 
farthest away from the origin are the ones which get the largest distortion relative to 
their physical counterparts. 

In real life, the number of available terms from the chiral expansion for the con- 
struction of a PA is very limited. Since the masses and decay constants of the first few 
vector and axial-vector resonances are known, one may envisage the construction of a 
rational approximant having some of its poles at the prescribed values given by the 
known masses of these resonances. If all the poles in the approximant are prescribed 
this way (as in the MHA), we have a PTA. On the contrary, when some of the poles are 
prescribed but some are also left free, then we have a PPA (see the previous section). 

Assuming that the first masses are known, let us proceed to constructing the PTAs 
(JIJJ). The lowest such PTA is T^iQ 2 ), which contains two poles at the physical masses 
of the p and the first A in the tower. Fixing the residue through the chiral expansion 
to be Co = —Fq, one obtains 

2 - F 2 M 2 p M\ 

Ts(g } = (Q 2 + M 2 )(Q 2 + M 2 ) ■ (24) 

Even though it has the same number of inputs (Co and the two masses), this rational 
approximant does not do such a good job as the PAs (fT5|) or (12TT) . For instance, C_ 4 
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is 2.3 times larger than the true value in Table [TJ As we have already stated, one 
way to intuitively understand this result is the following. The OPE is an expansion at 
Q 2 = oo and therefore knows about the whole spectrum because no resonance is heavy 
enough with respect to Q 2 to become negligible in the expansion, i.e., the infinite tower 
of resonances does not decouple in the OPE. Chopping an infinite set of poles down 
to a finite set may be a good approximation, but only at the expense of some changes. 
These changes amount to the appearance of poles and residues in the PA which the 
original function does not have. This is how the PA (fT9"j) manages to approximate 
the true function (JHJ). However, by construction, the PTA (|24|) does not allow the 
presence of any artificial pole because, unlike in a PA, all its poles are fixed at the 
physical values. Consequently, it only has its residues as a means to compensate for 
the infinite tower of poles present in the true function and, hence, does a poorer job 
than the PA (1191) . particularly in determining large-Q 2 observables like C_4. Indeed, 
the role played by the residues in the approximation can be appreciated by comparing 
the true values of the decay constants to those extracted from fl24l) . Although the one 
of the p is within 30% of the true value, that of the A is off by 100%. 

A different matter is the prediction of low-energy observables such as, e.g., the 
chiral coefficients. In this case heavy resonances make a small contribution and this 
means that the infinite tower of resonances does decouple!^! Truncating the infinite 
tower down to a finite set of poles is not such a severe simplification in this case, which 
helps understand why a PTA may do a good job predicting unknown chiral coefficients. 
Indeed, ( 1241) reproduces the value of C2 within an accuracy of 15%, growing to 22% in 
the case of C4. A global observable like I n averages the low and the high Q 2 behaviors 
and ends up differing from the true value (l2TJj) by 35%. This gives some confidence 
that observables which are integrals over Euclidean momentum may be reasonably 
estimated with MHA as, e.g., in the Bk calculation of Ref. [39J. 

Improving on the PTA (jUj) by adding in the first resonance mass from the vector 
tower produces the following approximant 

T i (n 2s_ a + bQ 2 J a = -13.5 x 10~ 3 , 

sW (Q 2 + M 2 )(Q 2 + M 2 A )(Q 2 + My) ' Wlth \ b = +1.33 x 10~ 4 , {Zb) 

where the values of the chiral coefficients Cq and C2 have been used to determine 
the parameters a and b. The prediction for C4 is much better now (only 2% off), 
in agreement with our previous comments. The prediction for C_4 is still very bad, 
becoming now 19 times smaller than the exact value. Nevertheless, it eventually gets 
much better if PTAs of very high order are constructed. For instance, we have found 
C_4 = —2.58 x 10 -3 for the approximant Tg with 9 poles. Similarly, we have also 
checked that the prediction of the chiral coefficients and the integral (TSUI) improve with 
higher-order PTAs. 

However, another matter is the prediction of the residues. For instance, the predic- 
tion for the decay constant of the state with mass My in ( 1251) is smaller than the exact 
value in the model (115p by a factor of 2. In general, we have seen that the residues 



14 This is because the residues F 2 in the Green's function (TT4|) stay constant as the masses grow. 
This behavior does not hold in the case of the scalar and pseudoscalar two-point functions [10] . 
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of the poles always deteriorate very quickly so that the residue corresponding to the 
pole which is at the greatest distance from the origin is nowhere near the exact value. 
We again explicitly checked this up to the approximant Tg, in which case the decay 
constant for this pole is almost 5 times smaller than the exact value. The conclusion, 
therefore, is that PTAs are able to approximate the exact function only at the expense 
of changing the residues of the poles from their physical values. Identifying residues 
with physical decay constants may be completely wrong in a PTA for the poles which 
are farthest away from the origin. 

As an intermediate approach between PAs and PTAs, there are the PPAs ([3]) where 
some poles are fixed at their physical values while some others are left free. The simplest 
of such rational approximants is i{Q 2 ) (see the previous section for notation). Fixing 
its 3 unknowns with M 2 , Cq and C 2 , one obtains 

P?,i(Q 2 ) = , n2 - ~ w* 2 - x , with r 2 R = 3.75 x lO" 3 , z R = 0.8665 . (26) 
(Q 2 + M 2 )(Q 2 + z R ) 

As can be seen, the mass (squared) of the first A resonance is predicted to be at z R 
which is sensibly smaller than the true value in ffT5lf^l . The rational function (1261) 
predicts C_4 = — r R = —3.75 x 10~ 3 which is a better determination than that of the 
PTA (jUj) with the same number of inputs, and C4 = —45.52 x which is not bad 
either. Concerning the pion mass difference, one gets 1^ = 5.22 x 10 -3 . However, as 
compared to the PAs fTiTJl) or (|2~T1) . the PPA f[2"6l) does not represent a clear improvement. 

In order to improve on accuracy of the PPA, one may try to use the mass and decay 
constant of the first resonance, M p and F p , in addition to the pion mass difference and 
the chiral coefficients Cq, C 2 and build the P 2 i{Q 2 ), which can be written as: 

F 2 M 2 a - F 2 M 2 Q 2 [ a = 17 43 x 1(T 3 

yWJ Q 2 + M 2 + (Q 2 + z c ) (Q 2 + z*) ' 1 z c = 1.24 + i 0.34 . {Z<) 

This PPA, upon expansion at large and small Q 2 , determines C_4 = —2.47 x 10^ 3 and 
C4 = —44.0 x 10~ 3 to be compared with the corresponding coefficient in Table [TJ The 
accuracy obtained is better than that of (I2"T|) . but this is probably to be expected since 
(|27|) has more inputs. 

Based on the previous numerical experiments done on the model in Eq. (I14|ll5p (and 
many others), we now summarize the following conclusions. Although, in principle, 
the PAs have the advantage of reaching the best precision by carefully adjusting the 
polynomial in the denominator to have some effective poles which simulate the infinite 
tower present in ffl~4l) . they have the disadvantage that some of the terms in the low- 
Q 2 expansion are required precisely to construct this denominator. This hampers the 
construction of high-order PAs and consequently limits the possible accuracy. 

When the locations of the first poles in the true function are known, there is the 
possibility to construct PTAs (with all the poles fixed at the true values) and PPA 
(with some of the poles fixed and some left free). As we have seen, although the PTA 
may approximate low-Q 2 properties of the true function reasonably well, the large-Q 2 



15 Intriguinly enough, this is also what happens in the real case of QCD [9| fTT] . 
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properties tend to be much worse, at least as long as they are not of unrealistically 
high order. The PPAs, on the other hand, interpolate smoothly between the PAs 
(only free poles) and the PTAs (no free pole). Depending on the case, one may choose 
one or several of these rational approximants. However, common to all the rational 
approximants constructed is the fact that the residues and/or poles which are farthest 
away from the origin are in general unrelated to their physical counterparts. 



4 The QCD case 

Let us now discuss the real case of large- N c QCD in the chiral limit. In contrast to the 
case of the previous model, any analysis in this case is limited by two obvious facts. 
First, any input value will have an error (from experiment and because of the chiral and 
large- A^ c limits), and this error will propagate through the rational approximant. And 
second, it is not possible to go to high orders in the construction of rational approxi- 
mants due to the rather sparse set of input data. In spite of these difficulties one may 
feel encouraged by the phenomenological fact that resonance saturation approximates 
meson physics rather well. 

The simplest PA to the function Q 2 U LR (— Q 2 ) with the right fall-off as Q~ A at large 
Q 2 is P 2 °(Q 2 ): 

^ = 1 + AQ> + BQ' ■ < 28 > 

The values of the three unknowns a, A and B may be fixed by requiring that this PA 
reproduces the correct values for F , L w and /„- [£| given by 

F = 0.086 ± 0.001 GeV , 
5m n = 4.5936 ± 0.0005 MeV I n = (5.480 ± 0.006) x 10" 3 GeV 4 , 

Li (0.5 GeV) < L w < L 10 (l.l GeV) L w = (-5.13 ± 0.6) x 10~ 3 . (29) 

The low-energy constant L w is related to the chiral coefficient C*2, in the notation of 
Eq. ([TBI) , by C 2 = — 4L 10 . Since L w does not run in the large- N c limit, it is not 
clear at what scale to evaluate L 10 (/x) [4"0] . In Eq. f[2"9"]) we have varied /i in the range 
0.5 GeV < fi < 1.1 GeV as a way to estimate 1/N C systematic effects. The central 
value corresponds to the result for L W (M P ) found in Ref. [12]. The other results in 
( 129]) are extracted from Refs. [H E] . 

Obviously, the PA (|28l) can also be rewritten as 

P °W 2 ) = 'JL ^ > ( 3 °) 

{Q 2 + z v ){Q 2 + z A ) 

in terms of two poles zy.A- In order to discuss the nature of these poles, we will define 
the dimensionless parameter ( by the combination 

C = -4L 10 -| = 2.06 ± 0.25 , (31) 



16 Recall that 1^ is, up to a constant, the electromagnetic pion mass difference Sm^ [17] and is 
defined in terms of IIlr as in Eq. ([20]) . 
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where the values in (I29I) above have been used in the last step. Imposing the constraints 
(12"9|) on the PA (150]) one finds two types of solutions depending on the value of (: for 
C > 2 the two poles z v ,a are real, whereas for ( < 2 the two poles are complex. At 
C — 2, the two solutions coincide. To see this, let us write the set of equations satisfied 
by the PA fl30]) as: 



y>2 



ZyZA 



-4L 10 = F£[— + — 

h = F log — . (32) 

z A - z v z v 

The first of these equations can be used to determine the value of the residue r 2 in 
terms of zyZA- In order to analyze the other two, let us first assume that both poles 
zy,A are real. In this case, they also have to be positive or else the integral I n will not 
exist because it runs over all positive values of Q 2 . Let us now make the change of 
variables 

z v = R (1 - x) , z A = R (1 + x) . (33) 

The condition z v ^a > translates into R > 0, \x\ < 1. In terms of these new variables, 
the second and third equations in (I3"2"|) can be combined into 

> 1, 1 + x . . 

C = - log , 34 

x 1 — X 

where the definition (|3T|) for ( has been used. With the help of the identity log(l + 
x/1 — x) =2 th _1 x (for |x| < 1), one can finally rewrite this expression as 

2 

( = — th _1 x , (x real) (35) 

which is an equation with a solution for x only if ( > 2. Once this value of x is found, 
the value of R can always be obtained from one of the last two equations ( 1321) and this 
determines the two real poles zy,A from (!33j) . 

On the other hand, when ( < 2, Eq. (|35|) does not have a solution. However, 
according to (13TT) . ( can also be smaller than 2. In order to study this case, we may 
use the identity th — 1 (i y) = i tan _1 (?/) to rewrite the above equation (135]) in terms of 
the variable x = i y (y real) as 

( = - tan" 1 ?/ , (y real). (36) 

y 

One now finds that this equation has a solution for y when C, < 2. In this case 
the poles of the PA (1281) are complex-conjugate to each other and can be obtained as 
zy t A — R(l±i y). These poles, obviously, cannot be associated with any resonance mass 
and this is why this solution has been discarded in all resonance saturation schemes 
up to now. However, from the point of view of the rational approximant (I28p there is 
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Co 


c 2 


c 4 


c 6 


c 8 


c_ 4 


C-6 


C-s 


-F 2 


-4Li 


-43 ± 13 


81 ±53 


-145 ± 120 


-4.1 ±0.5 


6±2 


-7 ±6 



Table 2: Values of the coefficients C<ik i n the high- and low-Q 2 expansions of 
Q 2 Tl LR (-Q 2 ) in Eq. (Q^ in units of 10~ 3 GeV 2 ~ 2k . Recall that C_ 2 = 0. 



nothing wrong with this complex solution, as the approximant is real and well behaved. 
From the lessons learned in the previous section with the model, there is no reason to 
discard this solution since, as we saw, rational approximants may use complex poles 
to produce accurate approximations. Therefore, we propose to use both the complex 
as well as the real solution for the poles zy,A, at least insofar as the value for ( ^ 2. In 
this case we obtain, using the values given in Eqs. ( 1291 . 

(C>2), r 2 = -(4.1 ±0.5) x 10~ 3 , z v = (0.77) 2 ± 0.15 , z A = (0.96) 2 ± 0.41 (37) 
(C < 2) , r 2 = -(3.9 ± 0.1) x 10~ 3 , z v = z* A = (0.66 ± 0.06) + % (0.25 ± 0.25) , (38) 

in units of GeV 6 for r 2 , and GeV 2 for zy,A- The two solutions in Eqs. (I37II38P have 
been separated for illustrative purposes only. It is clear that they are continuously 
connected through the boundary at ( = 2, at which value the two poles coincide and 
zv = z a — 0.72. The errors quoted are the result of scanning the spread of values in 
fl29|) through the equations (132]) . 

With both set of values in ( |37ll38l) . one can get to a prediction for the chiral and OPE 
coefficients by expansion in Q 2 and 1/Q 2 , respectively. These expansions of the PA can 
be done entirely in the Euclidean region Q 2 > 0, away from the position of the poles 
z v,Ai whether real or complex. Recalling the notation in Eq. (fl6l) . the above P^iQ 2 ) 
produces the coefficients for these expansions collected in Table [2J The values for the 
OPE coefficients C_4 i _ 6> _ 8 in this table are compatible with those of Ref. [17], after 
multiplying by a factor of two in order to agree with the normalization used by these 
authors. However, the spectrum in our case is different because of the complex solution 
in fl38|) . As we saw in the previous section with a model, this again shows that Euclidean 
properties of a given Green's function, such as the OPE and chiral expansions, or 
integrals over Q 2 > are safer to approximate with a rational approximant than 
Minkowskian quantities, such as resonance masses and decay constants. 

5 Conclusions 

In this article we pointed out that approximating large- iV c QCD with a finite number 
of resonances may be reinterpreted within the mathematical Theory of Pade Approxi- 
mants to meromorphic functions |18j . 

The main results of this theory may be summarized as follows. One may expect 
convergence of a sequence of Pade Approximants to a QCD Green's function in the 
large- N c limit in any compact region of the complex Q 2 plane except at most in a 
zero-area set [2D]. This set without convergence comprises the poles of the original 
Green's function together with some other artificial poles generated by the approximant 
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which the original function does not have. As the order of the PA grows, the previous 
convergence property implies that any given artificial pole either goes to infinity, away 
from the relevant region, or is almost compensated by a nearby zero. This symbiosis 
between a pole and a zero is called a defect. Although close to a pole the rational 
approximation breaks down, in a region which is far away from it the approximation 
should work well. 

We have reviewed the main results of this theory with the help of a model for the 
two-point Green's function (VV — AA). The simpler case of a Green's function of the 
Stieltjes type, such as the two-point correlator (VV), was previously considered in Ref. 
[2T] . We have seen in the case of this particular model how rational approximants 
create the expected artificial poles (and the corresponding residues) in the Minkowski 
region Re(g 2 ) > while, at the same time, yielding an accurate description of the 
Green's function in the Euclidean region Re(g 2 ) < 0. This happens in a hierarchical 
way: although the first poles/residues in a PA may be used to describe the physical 
masses/decay constants reasonably well, the last ones give only a very poor description. 
Therefore, it is in general unreliable to extract properties of individual mesons, such 
as masses and decay constants, from an approximation to large- N c QCD with only a 
finite number of states. Since a form factor, like a decay constant, is obtained as the 
residue of a Green's function at the corresponding pole(s), this also means that one 
may not extract a meson form factor from a rational approximant to a 3-point Green's 
function, in agreement with [16]. This observation may explain why the analysis of 
Ref. [43J, which is based on an extraction of matrix elements such as (tt\S\P) and 
(ir\P\S) from the 3-point function (SPP), finds values for the form factor which 
are different from those obtained in other analyses [S|. 

In spite of all the above problems related to the Minkowski region, our model 
shows how Pade Approximants may nevertheless be a useful tool in other regions of 
momentum space. We think that this is also true in the real case of QCD in the large- 
N c limit. In this case one may use the first few terms of the chiral and operator product 
expansions of a given Green's function to construct a Pade Approximant which should 
yield a reasonable description of this function in those regions of momentum space 
which are free of poles. In this construction, Pade Approximants containing complex 
poles, if they appear, should not be dismissed. 

We have also reanalyzed the simplest approximation to the (VV — AA) Green's 
function in real QCD which consists of keeping only two poles, and we have found 
that, depending on the value of the combination ( in Eq. fl3TT) . these two poles may 
actually be complex. 

However, if not all the residues and masses in a rational approximant are physical, 
this poses a challenge to any attempt to use a Lagrangian with a finite number of 
resonances such as, for example, the ones in Ref. [21 E], for describing Green's functions 
in the large--/V c limit of QCD. Even if these Lagrangians are interpreted in terms of 
PTAs, with the poles fixed at the physical value of the meson masses, we have seen how 
the residues then get very large corrections with respect to their physical counterparts. 
Can these residues be efficiently incorporated in a Lagrangian framework? We hope to 
be able to devote some work to answering this and related questions in the future. 
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APPENDIX 



Here we will show how the PAs constructed from the OPE do not in general re- 
produce even the first resonances in the spectrum, unlike those constructed from the 
chiral expansion. Again, we will use the model of section 3 as an example. Recalling 
the definition of the OPE given in Eq. ffl6l) . with the corresponding coefficients (ITT)) , it 
is straightforward to construct a PA in 1/Q 2 around infinity, i.e. by matching powers 
of the OPE in 1/Q 2 . The construction parallels that in Eq. fl2j) but with the replace- 
ment z = 1/Q 2 - Since the function Q 2 I1lr{— Q 2 ) behaves like a constant for Q 2 — > 0, 
we will consider diagonal Pade Approximants, i.e. of the form P"(l/Q 2 ), in order to 
reproduce this behavior. Figure [2] shows the position of the poles and zeros of the PA 
P|q (— 1/q 2 ) in the complex q 2 plane. As it is clear from this plot, the positions of the 
poles have nothing to do with the physical masses in the model, given by Eqs. (II lti!5j) . 
even for the lightest states. This is to be contrasted with what happens with the PA 
constructed from the chiral expansion around Q 2 , which is shown in Fig. [TJ The differ- 
ence between the two behaviors is due to the fact that, while the chiral expansion has 
a finite radius of convergence, the radius of convergence of the OPE vanishes because 
this expansion is asymptotic. 
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